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Approach-to-equilibrium molecular dynamics simulations have been used to study thermal trans¬ 
port in nanocrystalline graphene sheets. Nanostructured graphene has been created using an it¬ 
erative process for grain growth from initial seeds with random crystallographic orientations. The 
resulting cells have been characterized by the grain size distribution based on the radius of gyration, 
by the number of atoms in each grain and by the number of atoms in the grain boundary. Intro¬ 
duction of nanograins with a radius of gyration of 1 nm has led to a significant reduction in the 
thermal conductivity to 3% of the value in single crystalline graphene. Analysis of the vibrational 
density of states has revealed a general reduction of the vibrational intensities and broadening of 
the peaks when nanograins are introduced which can be attributed to phonon scattering in the 
boundary layer. The thermal conductivity has been evaluated as a function of the grain size with 
increasing size up to 14 nm and it has been shown to follow an inverse rational function. The grain 
size dependent thermal conductivity could be approximated well by a function where transport is 
described by a connection in series of conducting elements and resistances (at boundaries). 


I. INTRODUCTION 

Graphene is a material of great interest from a tech¬ 
nological point of view owing to its high carrier mobility 
and optical transparency. Furthermore, it provides ex¬ 
cellent physical strength and is chemically inert towards 
many environmentally ubiquitous speciesSuch prop¬ 
erties paved the way for graphene’s career in applica¬ 
tions such as field effective transistors, Schottky junction 
diodes or as transparent and flexible displays in various 
electronic devices 

Miniaturization of electronic devices is an important 
topic in recent research. A major issue for such devices 
is the required heat dissipation which becomes more and 
more important with decreasing dimensions. Crystalline 
graphene is a promising material in this regard. It inheres 
both high electronic carrier mobility and excellent ther¬ 
mal conductivity which additionally provides the func¬ 
tion as heat dissipator. 

Fabrication of single-crystalline graphene sheets, how¬ 
ever, is not trivial. Common fabrication methods in¬ 
clude epitaxial film growth and chemical vapor deposi¬ 
tion (CVD), where the latter is promising in particular 
for large-scale production. In all methods, defects are in¬ 
troduced in the crystalline structure resulting from lim¬ 
itations of the kinetic properties in the growth process 
and defects of the substrates. These defects limit the 
size of single-crystalline domains rendering the graphene 
sheets rather polycrystalline. The introduction of grain 
boundaries and defects in such polycrystalline graphene 
can drastically change electronic and thermal transport 
properties of the system. In fact, previous experimental 
measurements of the thermal conductivity in graphene 
resulted in values ranging from 600 to 10000 W/mK.^ 
This wide range in the value of thermal conductivity 
demonstrates its high sensitivity to various process con¬ 
ditions, sample size, and structural features including de¬ 


fects and grain boundaries. 

Characterization and control of the structural proper¬ 
ties of polycrystalline graphene, such as the grain size, 
and their effect on the electrical and thermal trans¬ 
port properties are thus of major interest in research of 
graphene in electronic and thermoelectric devices. Re¬ 
garding the effect of environmental parameters, a sig¬ 
nificant difference in the thermal conductivity has been 
found for suspended (^1600 W/mK)pIl and Si 02 sup¬ 
ported (600 W/mK)P^ graphene. Furthermore, thermal 
conductivity of CVD produced graphene has been shown 
to be reduced by 30% when wrinkles are introduced in 
the structure.^ A systematic investigation of the size 
effect of single-crystalline domains on the thermal con¬ 
ductivity of graphene based on experiments, however, is 
still matter of investigation. 

Due to the complexity of generating trustworthy poly¬ 
crystalline simulation cells, little research has been done 
so far on the theoretical investigation of the properties 
of polycrystalline graphene. Applying atomistic simu¬ 
lations it has been shown recently that the mechani¬ 
cal strength of graphene is reduced by 50% whengrain 
boundaries in form of polycrystals are introduced.^ Fur¬ 
thermore, equilibrium molecular dynamics (EMD) simu¬ 
lations have been carried out where the thermal conduc¬ 
tivity has been described as a function of the polycrys¬ 
talline grain sizes 

In this work, in addition to the grain size dependent 
thermal conductivity of nanocrystalline graphene, the 
morphology of crystalline domains has been analyzed in 
detail. Furthermore, the phonon spectrum has been ex¬ 
tensively investigated based on the vibrational density 
of states (VDOS) and the accumulation of bulk ther¬ 
mal conductivity. First, we developed a model to create 
nanocrystalline graphene sheets with different grain sizes 
and determined their particle size distributions. Ther¬ 
mal transport properties of these systems have then been 
studied by approach to equilibrium molecular dynamics 
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(AEMD) calculations Two specific systems with a 
grain size of 1 and 2.5 nm, respectively, have been cho¬ 
sen for a detailed elaboration of the thermal transport 
including the extrapolation of the bulk thermal conduc¬ 
tivity and the analysis of the phonon spectrum based on 
the accumulation function of the thermal conductivity. 
Furthermore, the vibrational density of states (VDOS) 
in nanocrsytalline graphene has been evaluated and com¬ 
pared to the VDOS of pristine graphene. Finally, the 
thermal conductivity of systems with varying average 
grain size from 0.7 to 14 nm has been investigated show¬ 
ing increasing thermal conductivity with increasing grain 
size. The thermal conductivity has been approximated as 
rational functions both of the grain size and the number 
(concentration) of atoms in the grain boundaries. 


II. METHODS 



A. Generation of nanocrystalline graphene 


In this study, thermal transport in extended (bulk) 
nanocrystalline graphene has been studied using atom¬ 
istic simulations. Computer generation of nanostruc- 
tured materials, however, is not trivial, in particular, 
when periodic boundary conditions have to be met as it is 
the case in extended bulk systems. In analogy to exper¬ 
iments, production of such nanostructured systems can 
be realized by bottom-up or top-down methods. Top- 
down methods include the creation of a homogeneous 
amorphous or crystalline simulation cell which is used 
as a matrix material. Manual manipulation by scissor¬ 
ing and insertion of different seeds and structures is then 
used to create nanostructured materials. This method 
has been used, for example, for the creation of nanocrys¬ 
talline SiGe alloys.^ In contrast, in bottom-up meth¬ 
ods nanostructures are constructed from single atoms 
or small building blocks (for example single-crystalline 
domains). Two-dimensional (2D) graphene sheets simu¬ 
lated here have been simulated using the bottom-up ap¬ 
proach. 

An algorithm has been developed where, for a given 
rectangular simulation cell with width and length 
(orthogonal to and in direction of the thermal transport, 
respectively), a specified number of differently oriented 
crystalline grains has been placed inside the simulation 
box (Figure Q. Initial nucleation sites have been placed 
randomly in the simulation cell with the constraint that 
neighboring seed atoms have a distance dmin between 
each other of at least one fifth of the ideal average grain 
size «i„ = %- = lyilgl). 

An arbitrary angle between 0 ° and 60 ° has been as¬ 
signed to each nucleation site to define the orientation of 
the crystalline growth of graphene. Under-coordinated 
atoms on the edge of each grain have been added to a 
list of possible reactive sites. In an iterative process, the 
reactive site to which the next atom has been added, has 


FIG. I. Schematic of the simulation cell and the algorithm 
developed for creation of nanocrystalline graphene. In the 
above scheme, cn{Cj^i) indicates the coordination number of 
the respective atom i. 


been randomly selected from this list. The list of reactive 
atoms has been updated in each step until no reactive 
atoms were left. Grain growth has been terminated, i 
e. atoms have been eliminated from the list of reactive 
atoms, when the respective C atom was three-fold coor¬ 
dinated or atoms from neighboring cells were closer than 
95% of the equilibrium C-C bond length (here: 1.3978 

A). 


Creation of simulation cells by this method only con¬ 
siders geometrical aspects. Interatomic interactions are 
neglected. Furthermore, simulation cells are created only 
in two dimensions, omitting corrugation in off-plane di¬ 
rection. To account for the latter, simulation cells have 
been relaxed during molecular dynamics (MD) simula¬ 
tions at high temperatures. Details on the MD simula¬ 
tions can be found in Section IIB Velocity rescaling has 
been applied for this equilibration process using a time 
step of 0.5 fs. It has been initiated at a temperature of 
300 K for 10 ps. Subsequently, the cells have been heated 
to a temperature of 4000 K with a heating rate of 105 
K/ps (in 35 ps) and relaxed at this temperature for 255 
ps. The cells have been cooled back to 300 K with a 
cooling rate of 10 K/ps (in 370 ps). This procedure has 
been followed by a microcanonical (NVE) run for 200 ps. 


Several geometrical criteria have been used to charac¬ 
terize the grain size of the created nanocrystalline sam¬ 
ples. One possible way is the determination of an ideal 
radius rid assuming the particles to have a circular shape 
calculated from the area of the grain which is estimated 


from the number of atoms in the grain (rid 



However, this characterization does not consider differ¬ 
ent shapes of the grains. Therefore, we also included the 
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FIG. 2. Simulation cells of nanocrystalline graphene with a 
mean radius of gyration of (a) 5.1, (b) 2.5 and (c) 1.4 nm. 
Atoms assigned to the grain boundary are indicated in blue. 


characterization by the radius of gyration which 

has been determined for each seed j according to Eq. 
where rij is the position of atom i in grain j and rcMj 
and nc,j are the center of mass and the number of atoms 
in the grain, respectively. 
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At least two different configurations have been created 
for each sample dimension and grain size. For systems 
with an average radius of gyration of 1 nm, the trans¬ 
verse section orthogonal to the direction of the heat 
flow has been changed from 5 to 20 nm. Furthermore, 
for such systems, the sample length Lz in direction of the 
thermal transport has been changed from 50 to 1000 nm. 
The effect of the sample length has been investigated ad¬ 
ditionally for samples with rG=2.5 nm and Lx=20 nm. 
These calculations have been carried out to verify the re¬ 
liability of the used simulation cells for the determination 
of the thermal conductivity. Based on these results, the 
average radius of gyration has been changed from 0.7 nm 
to 14 nm using samples with a cell length of 200 nm. 

Atoms belonging to the grain boundaries have been 
identified based on their energy (energy/atom). For this 
purpose, the energy/atom distribution has been evalu¬ 
ated for initial configurations of each sample. It shows a 
main peak for the atoms belonging to the center of the 
grains with a shoulder evolving at higher energies corre¬ 
sponding to atoms in the grain boundaries. Atoms with 
energies higher than 2% of the main peak value have been 
assigned to the grain boundaries. By this technique some 
atoms from the center of the grain might be selected as a 
result of corrugation in out-of-plane direction. However, 
the amount is negligible. Nanocrystalline graphene grain 
boundary atoms selected in this manner are shown in 
Figure for samples with grain sizes of rG=l-4, 2.5 and 
5.1 nm. This characterization allowed us to determine 
the concentration of grain boundary atoms cgb in each 
sample. 


B. AEMD simulations 

Molecular dynamics si mulat ions have been performed 
using the lammps code.^^^^ Covalent interactions be¬ 


tween carbon atoms have been described by the second- 
generation reactive empirical bond order (REBO) poten¬ 
tial which has previously been sho wn to give reasonable 
results for graphene-based systems.l^^^^ 

The thermal conductivity has been determined based 
on the AEMD methodology. In this formalism, the sim¬ 
ulation cell is firstly divided into two regions with equal 
length Lzl2. One of these two compartments is equi¬ 
librated at a high temperature (Th=400 K), the other 
compartment at a low temperature (Tc=200 K) by veloc¬ 
ity rescaling. This creates an initial step-like function of 
the temperature along the sample length in z-direct ion. 
Next, the evolution of the average temperature in the hot 
(Th) and cold (Tc) reservoir is recorded during a transient 
regime towards equilibrium of microcanonical evolution. 
Based on Fourier’s theorem of thermal transport and the 
given step-like initial temperature profile, the evolution 
of the temperature gradient (AT = Th — Tc) follows 
AT (t) = where R = ^ is the thermal 

diffusivity with the density p of the material and its spe¬ 
cific heat Cy. This expression is fitted to the temperature 
gradient obtained from the simulations to determine ti. 
More det ails o n the AEMD methodology can be found 
elsewhere 


III. RESULTS AND DISCUSSION 

A. Morphology of nanocrystalline graphene 
samples 

Nanocrystalline graphene samples have been charac¬ 
terized by the number of C atoms nc in each grain and 
by their radius of gyration vq which has been determined 
according to Eq. The particle size distribution (PSD) 
based on nc is best described by a lognormal distribu¬ 
tion as shown in Figure on the example of samples with 
average vq of 1 and 2.5 nm. The data has been fitted to 
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resulting in a mean number of atoms per grain = 
e^’^) of 237 and 1419, respectively. The PSD based on 
the radius of gyration, however, is better described by a 
normal distribution and accordingly has been fitted to 
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For a sufficiently large number of grains, the mean ra¬ 
dius of gyration obtained from Eq. is equivalent 
to its arithmetic mean. Small variations, however, have 
been observed when the number of seeds is respectively 
low and not sufficient for a statistical analysis. In a sim¬ 
ulation cell with dimensions 40 nm x 200 nm and a total 
number of 200 grains {re = 2.5 nm, Figure]^, distributed 
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FIG. 3. PSD of nanocrystalline graphene with an average 
radius of gyration re of (a,b) 1 and (c,d) 2.5 nm. The PSD 
follows a logarithmic normal distribution as a function of the 
number of atoms per grain (a,c) and a normal distribution as 
a function of re (b,d). 



FIG. 4. Normal distribution of nanocrystalline samples with 
average radii of gyration from 0.7 to 10 nm. 


in 50 bins, the mean of the normal distribution is slightly 
lower (2.38 nm) than the arithmetic mean. The fitted 
normal distribution of simulated samples with Lz=2{){) 
nm and varied grain sizes from 0.7 to 10 nm are shown 
in Figure]^ 

Grain boundaries have been further analyzed based on 
the number (concentration) of atoms in the grain bound¬ 
ary cgb- Assuming ideal circular shape of the grains, cqb 
is related to the particle radius rid according to 


Cgb 


Agb _ ^ t 
Ag Gd 


( 4 ) 


where Aq is the total the area of one grain and Aqb is 
the area of the grain boundary with a thickness t. Fitting 
Cgb to eq. results in a thickness of 0.164 nm (Fig¬ 
ure]^. From the density (atoms/area) of various samples 
calculate here we estimated the spacial extension of one 
atomic layer to be 0.16 nm. We therefore conclude that 
the grain boundaries in nanocrystalline graphene consist 
in average of one atomic layer of C atoms. 
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FIG. 5. Number of atoms in the boundary region as a function 
of the radius rid, assuming ideal circular shape. 



B. Sample size effects 

1. Length-dependent thermal eonduetivity 

The effect of the sample length Lz on the estimated 
thermal conductivity has been investigated for samples 
with average grain size tq of 1 and 2.5 nm. The sam¬ 
ple length has been increased from 50 nm to 1 qm. The 
size effect on thermal conductivity of crystalline graphene 
is under major dispute in recent literature. Several ex¬ 
perimental and theoretical studies have claimed diver¬ 
gence of the thermal conductivity with the length of the 
sample.l^^M^ However, more recently it has been proven 
that, while thermal conductivity grows logarithmically 
with the sample length up to ^10 pm, it is definitely up¬ 
per limited in samples long enough to allow a diffusion 
transport regime f or both single and collective phonon 
excitations 

Accordingly, here, we describe the length dependent 
thermal conductivity of crystalline and nanocrystalline 
graphene with the linear reciprocal (rational) function 


^oo 



from which the convergent bulk thermal conductivity 
at infinite sample size can be estimated. 

In the studied range of sample sizes both in crystalline 
and the two nano cyst alline samples, the 1//^ to l/L^ be¬ 
havior follows a linear trend according to Eq.|^(Figure]^ 
indicating that phonon properties in such materials can 
be approximated well by an average value.l^ A bulk ther¬ 
mal conductivity hZoo of 60.2 and 26.6 W/mK (Table 
has been estimated for nanocrystalline graphene with vq 
of 2.5 and 1 nm, respectively. These values correspond 
to 7 and 3%, respectively, of the bulk thermal conductiv¬ 
ity of crystalline graphene (859.7 W/mK). The relative 
change of the thermal conductivity with respect to the 
one of pristine graphene is in excellent agreement with a 
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TABLE L Optimized parameters Koo and A of the approxima¬ 
tion function Eq. [^for the thermal conductivity of nanocrys¬ 
talline graphene with average grain size of 1 and 2.5 nm and 
of pristine graphene as a function of the sample size. 


re [nm] 


A [nm] 

Lx [nm] 

1 

25.2±0.3 

27.7±3.0 

5 

1 

23.9±0.1 

22.4±1.1 

10 

1 

26.6±0.6 

38.3±5.8 

5-80 

2.5 

60.2±0.7 

40.4±3.3 

5-80 

crystalline 859.7±29.6 455.9±40.0 

3-10 


recent study where the thermal conductivity of polycrys¬ 
talline graphene has been estimated using equilibrium 
molecular dynamics simulations.!!^ 



1/Lqnm'] 

EIG. 6. Inverse of thermal conductivity 1/k, as a function 
of the sample length l/L^ of nanocrystalline graphene with 
average grain size re = 1 and 2.5 nm. 


2. Effect of transverse section 

The effect of the transverse of the samples on the esti¬ 
mated thermal conductivity has been evaluated for sam¬ 
ples with average grain size of rQ=l nm where an in¬ 
tersection of 5 and 10 nm has been simulated. Several 
configurations have been considered at each length. 

For an intersection of 5 nm, results in 25.2 W/mK 
and the characteristic length A is estimated with 27.7 
nm. Increase of the intersection had only little effect on 
the thermal conductivity. At an intersection of 10 nm, 
/^oo and A have been determined to be 23.9 W/mK and 
22.4 nm, respectively (Table [I|. We thus consider an in¬ 
tersection of 5 nm to be sufficiently large for a reasonable 
estimation of n. 


C. Accumulation of thermal conductivity 

The spectral contribution of phonons to the thermal 
conductivity has been estimated by an accumulation 
function. The accumulation function describes the ratio 
of the estimated thermal conductivity at a certai n sam¬ 
ple length to the bulk thermal conductivity.!^!!^ It can 
thus be understood as the contribution to the thermal 
conductivity of phonons that have a MFP (A) smaller 
than the corresponding sample length Lz according to 



where / (A) is the distribution function of the MFP of 
dominating phonons. 

A normal distribution has been assumed for /(A), 
where A enters into the equation as the logarithm of Lz 
and has been normalized to the characteristic length A 
(A = log ("x))- Following Ref. we set 


The characteristic length A, on the other hand, changes 
only marginally with the grain size. It has been estimated 
to be 40.4 and 38.3 nm for grains with tq of 2.5 and 1 nm, 
respectively. This indicates similar convergence behavior 
of the thermal conductivity in the two samples. In fact, 
according to Eq. 90% of the bulk thermal conductivity 
are obtained at a sample length Lz of 364 nm and 345 
nm for a grain size of 2.5 nm and 1 nm, respectively. In 
contrast, in crystalline graphene, samples with a length 
of more than 4 qm (A = 455.9 nm) have to be simulated 
to reach 90% of the bulk thermal conductivity. 

These results demonstrate that the mean free path 
(MFP) of dominant phonons present in crystalline 
graphene is drastically reduced when grain boundaries 
in the form of nanosized grains are introduced. The size 
of the grains, however, influences the MFP of dominant 
phonons only marginally, at least for the grain sizes stud¬ 
ied here (2.5 and 1 nm). 


/(A) = 


1 




exp 




( 7 ) 


The parameters fi and a of the normal distribution 
of the phonon MFP have been estimated fitting the ac¬ 
cumulation of thermal conductivity as a function of the 
sample lengths Lz to Eq. The statistical value jl can 
be translated into //mfp = A • 10^ which corresponds to 
the average value of the phononic MFP of the system. 


k{A) 


1 
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1 + erf 
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Figure shows the accumulation function that has 
been optimized for both nanocrystalline systems in com¬ 
parison to the results of crystalline graphene. The corre¬ 
sponding normal distribution is represented in Figure 
along with the optimized parameters in Table [llj 
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FIG. 7. (a) Accumulation and (b) distribution function of 
nanocrystalline and crystalline graphene. 


TABLE II. Estimated values of variance and average mean 
free path /xmfp of dominant phonons. 


re [nm] 

//MFP [nm] 

0-2 

1 

30.46 ± 0.107 

0.718 0.978 

2.5 

33.46 ± 0.036 

0.686 0.997 

crystalline 451.38 ± 0.015 0.535 0.995 


The average MFP of dominant phonons //mfp in 
nanocrystalline graphene with a grain size of 2.5 and 1 
nm is estimated to be 30.5 and 33.5 nm, respectively, in¬ 
dicating that decrease of the grain size from 2.5 to 1 nm 
hardly changes the MFP spectrum of dominant phonons. 
In both cases, a significant effect of the grain size on the 
bulk thermal conductivity is observed. We therefore ar¬ 
gue that lattice vibrations in the nanocrystalline mate¬ 
rial are similar to the vibrations in pristine graphene but 
with additional contributions from phonon scattering at 
the grain boundaries. 


D. Analysis of the vibrational density of states 

In order to better characterize the phonon spectrum 
in nanocrystalline graphene, the vibrational density of 



FIG. 8. Vibrational density of states (VDOS) of crys¬ 
talline and nanocrystalline graphene in the range of acoustic 
phonons. 


states (VDOS) has been estimated for several sam¬ 
ples including crystalline graphene and nanocrystalline 
graphene with rG=0.7, 1, 2.5 and 5.1 nm. Furthermore, 
the VDOS has been calculated specifically of atoms in the 
grain boundary. The VDOS has been calculated as the 
Fourier transform of the velocity autocorrelation function 
where the atomic velocities have been collected during a 
canonical simulation (NVT) at 300 K for 10 ps preceded 
by a canonical equilibration at 300 K for 200 ps. The ac¬ 
curacy of this method has been compared to results ob¬ 
tained from density functional theory (DFT) calculations 
through the diagonalization of the dynamical matri:jP^ 
and have been shown to be in general agreement. This 
is in particular the case in the acoustic range of phonon 
modes (<50THz) which are mainly contributing to ther¬ 
mal transport in graphene verifying the reliability of our 
model to describe the VDOS in graphene-based materi¬ 
als. 

Figure shows the VDOS of acoustic phonons of pris¬ 
tine and nanocrystalline graphene and of the atoms be¬ 
longing to the grain boundary. The main peaks of acous¬ 
tic phonons in crystalline graphene are observed at ^12 
and ^22 THz. The peak at 22 THz is extended by a 
shoulder to higher frequencies resulting from an addi¬ 
tional peak at ^28 THz. In nanocrystalline graphene, 
an overall broadening of all peaks is observed combined 
with a decrease in intensity. In particular, a shoulder 
at lower frequencies of the 12 THz peak seems to evolve 
which indicates enhanced vibrations at ^9 THz. Simi¬ 
larly, the shoulder at higher energies of the 22 THz peak 
which is most pronounced at ^28 THz flattens out in 
nanocrystalline samples which can possibly result from 
additional vibrations at ^25 THz. These observations 
are confirmed regarding the VDOS of grain boundary 
atoms (black line). There, the main vibrations of crys- 
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talline graphene at 12 and 22 THz disappeared due to the 
absence of any ordered crystalline structure. Instead, the 
density of vibrations increases monotonously from 0 to 22 
THz and is basically unchanged for frequencies from 22 
to 30 THz. This indicates enhanced vibrational states 
between 15 and 20 THz and between 25 and 28 THz 
with respect to crystalline graphene. In summary, we 
consider the decrease in intensity and the broadening of 
the peaks as the main reasons for the observed decreased 
in the thermal conductivity of nanocrystalline graphene. 
Whereas phonon scattering at the grain boundaries can 
be assigned to vibrations between 15 and 20 THz and 
between 25 and 28 THz. 


E. Grain size effects 


In the interest of computational costs, the effect of the 
grain size in nanocrystalline graphene on its thermal con¬ 
ductivity has been studied at a finite sample length of 
I/; 2=200 nm. 

For an intersection of 20 nm, samples with initial seed 
densities from 0.025 to 0.32 seeds/nm^, resulting in grain 
sizes re from 0.7 to 2.5 nm, have been simulated. To 
determine the thermal conductivity of larger grains, the 
intersection had to be increased to 40 nm and 80 nm. Ini¬ 
tial seed densities have been reduced down to 7.5x10“^ 
seeds/nm^ leading grains sizes up to 14 nm. Following 
the results of the previous paragraphs, marginal devi¬ 
ation in the estimated thermal conductivity is expected 
when the intersectional length is varied. Thus, results 
of the thermal conductivity of samples with different 
will be compared directly. 

Thermal conductivity of nanocrystalline graphene as a 
function of the grain size is shown in Figurej^ It increases 
continuously from 18.3 to 156 W/mK for increasing grain 
size (re) from 0.7 to 14 nm. The introduction of grain 
boundaries in the crystalline graphene can be regarded 
as a connection in series of resistances. Based on this a 
rational function is used to describe the dependency of n 
on re (Eq.[^ where /^c-Graph is the thermal conductivity 
of crystalline grains which is assumed to be constant. 
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A similar argumentation has been used previ¬ 
ously where the thermal conductivity of polycrystalline 
graphene has been calculated using equilibrium molecu¬ 
lar dynamics.^ 

The parameters /^c-Graph and Rqb of Eq. [^have been 
optimized with 282 W/mK and 8.2x10“^^ rrrK/W, re¬ 
spectively. It has to be noted that the value of the ther¬ 
mal resistance might not be representative since it is very 
sensitive to the value of vq which is only a rough esti¬ 
mation of geometrical properties of the grains. It gives 


Tq [nm] 

0 2 4 6 8 10 12 14 16 
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FIG. 9. Thermal conductivity of nanocrystalline graphene 
as a function of the concentration of atoms in the grain 
boundary cqb (solid line) and of the radius of gyration re 
(dotted line). 


indications for the thermal resistance but should rather 
be regarded as a fitting parameter for the description of 
K 8iS 8i function of vq. The constructed constant thermal 
conductivity /^c-Graph slightly overestimates the thermal 
conductivity of pristine graphene at this cell length (262 
W/mK) by 7.5 %. /^c-Graph thus overestimates the ac¬ 
tual value of crystalline graphene at this cell length by 
24%. 

This is striking since it predicts that the thermal 
conductivity of the crystalline domains (with a length 
smaller than the simulation cell) would be higher than 
the reference value of the thermal conductivity of crys¬ 
talline graphene in a simulation cell with the same di¬ 
mensions. However, it can partly be explained by the 
fundamentals of the approximation function of hi{rG)- 
In Eq. the number of grain boundaries contributing to 
the thermal resistance Rqb is assumed to be the same as 
the number of grains riseeds? instead it should be reduced 
by one {ugb = ^seeds ~ !)• This approximation is valid 
for a reasonably large number of grains ^ 1). 

However, since simulations have been carried out at a 
finite sample length Lz, '^seeds decreases with increasing 
tq. Thus Eq. [^overestimates the actual thermal conduc¬ 
tivity for large re¬ 
in Eq. [^ it is assumed that crystalline domains are 
spherical with a radius vq and the length of the bound¬ 
aries perpendicular to thermal transport, effectively con¬ 
tributing to the thermal resistance Rqb^ is equal to the 
diameter (2rGB) of the grains. To account for the actual 
geometry of grains, the thermal conductivity has been 
described as a function of the number (concentration) of 
atoms in the grain boundaries (cgb) according to 


(cgb) 


1 + 
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where di corresponds to the thickness of one atomic 
layer of C atoms. Eor the samples calculated here an 
average value of 0.16 nm has been used. Applying this 













expression for the overall thermal conductivity n of the 
nanocrystalline simulation cell, the thermal boundary re¬ 
sistance i?GB is approximated with 2.8 x 10-11 m^K/W 
and the constant thermal conductivity of the crystalline 
domains /^c-Graph is found to be 269 W/mK. As men¬ 
tioned before, i?GB will be regarded only as a fitting 
parameter rather than the actual physical value of the 
thermal resistance. The constant thermal conductivity 
^c-Graph, oii the other hand, can directly be compared 
to the thermal conductivity of crystalline graphene. It 
still overestimates the value of crystalline graphene (262 
W/mK), however, in this case only by 2.6%. This demon¬ 
strates the relevance to account for the actual properties 
of the grain boundaries (density of such) when estimating 
the thermal transport. 

Both approximations indicate that the thermal con¬ 
ductivity of the crystalline domain /^c-Graph is estimated 
well by its bulk-like value. It demonstrates that the ther¬ 
mal conductivity of nanostructured systems such as poly- 
crystalline graphene can be estimated considering grain 
boundaries and crystalline parts as a connection of re¬ 
sistances in series. In this scheme phonon scattering ef¬ 
fects are included in the thermal resistance Rqb of the 
grain boundaries and the thermal conductivity of crys¬ 
talline domains can be described by the thermal conduc¬ 
tivity of the crystal without grain boundaries. However, 
here, grain boundaries are relatively small and can be 
described as infinitely small. It is supposed that for sys¬ 
tems with larger grain boundaries, phonon scattering ef¬ 
fects cannot only be included in the thermal boundary 
resistance but need to be included as a separate length 
dependent resistance/conductivity. 

IV. CONCLUSIONS 

An iterative algorithm has been implemented to create 
samples of 2D nanocrystalline graphene sheets respecting 
periodic boundary conditions. The average grain size of 
single-crystalline domains has been controlled by the ini¬ 
tial density of seeds per unit area. The grain size has 
been characterized by the radius of gyration, the num¬ 
ber of particles per atom and the concentration of atoms 
belonging to the grain boundaries. The grain size dis¬ 
tribution of the nanocrystalline materials described by 
the radius of gyration follows normal distribution. Char¬ 
acterizing the size of particles by their area, i.e. num¬ 
ber of particles per grain, a lognormal distribution of 
the particle size has been found. The thickness of grain 
boundaries has been estimated from the change in con¬ 
centration of grain boundary atoms as a function of the 
idealized circular radius. It has been estimated to be 0.16 
nm corresponding to one atomic layer. 

Thermal transport of such nanocrystalline graphene 
sheets has been exploited using AEMD simulations at an 
average temperature of 300 K. The sample size depen¬ 


dence of the thermal conductivity has been investigated 
in nanocrystalline samples with an average grain size of I 
and 2.5 nm showing a linear trend between 1/k and I/T^. 
Based on this, the bulk thermal conductivity at infinite 
sample length could be estimated. In a sample with aver¬ 
age grain size of I nm the bulk thermal conductivity has 
been estimated to be 26.6 W/mK which corresponds to 
3% of the thermal conductivity of crystalline graphene. 

Based on the accumulation function of the thermal 
conductivity, the average mean free path of phonons con¬ 
tributing to thermal conductivity has been estimated for 
both nanocrystalline samples and crystalline graphene. 
Introduction of grain boundaries resulted in a remark¬ 
able decrease of the estimated mean free path of crys¬ 
talline graphene (451 nm) to 30 nm. The grain size, 
however, had only marginal effect on the average mean 
free path. 

Analysis of the vibrational density of states (VDOS) of 
crystalline and nanocrystalline graphene showed a gen¬ 
eral decrease in the density of vibrational states and a 
broadening of the peaks when grain boundaries in form 
of nanograins are introduced. The VDOS of grain bound¬ 
ary atoms in particular revealed increased vibrations be¬ 
tween 15 and 20 Hz and between 25 and 28 Hz which 
can be assigned to phonon scattering in the boundary 
layer and lead to the observed broadening of the VDOS 
in nanocrystalline graphene. 

Thermal conductivity at a sample length of 200 nm 
has been analysed as a function of the grain size. In 
the investigated range {vq = 0.7 to 14 nm), the thermal 
conductivity can be approximated well with a rational 
function based on the idea of grain boundaries to behave 
as thermal resistances connected in series. 

A better description of the overall thermal conductiv¬ 
ity of the nanostructured system, however, is obtained 
when the concentration of grain boundary atoms is con¬ 
sidered. It is proposed that for such nanostructured sys¬ 
tems, the overall thermal resistance can be described well 
as a connection of series where phonon scattering effects 
are included solely in the resistance of the grain bound¬ 
ary and the conductivity of the crystalline domains can 
be estimated by the conductivity of a crystal without 
grain boundaries. 
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